Mon. Not. R. Astron. Soc. 0Q0.[TinT](2Q08) Printed 23 September 2008 (MN KTeX style file v2.2) 



Integral Field Spectroscopy of HH 262: The Spectral Atlas 



oo 
o 
o 

m 

6 



> 

\o 
a^ 
rn 

o 
oo 
o 



R. Lopez/* B. Garcia-Lorenzo ^* S. F. Sanchez,^* G. Gomez,^'^* R. Estalella/* and 
A. Riera^* f 

^ Departament d'Astronomia i Meteorologia, Universitat de Barcelona, Marti i Franques 1, E-08028 Barcelona, Spain 
^Instituto de Astrofisica de Canarias, E-38200 La Laguna, Spain 

^Centro Astronomico Hispano-Alemdn de Calar Alto, Jesus Durban Remon 2-2, E-04004 Almeri'a, Spain. 
■^GTC Project Office, GRANTECAN S.A. (CALP). E-38712 Breiki Baja, LaPalma, Spain. 

^Dept. Fi'sica i Enginyeria Nuclear. EUETI de Barcelona. Universitat PoUtecnica de Catalunya. Compte d'Urgell 187, E-08036 Barcelona, Spain 



X 
S3 



Accepted . Received ; 



ABSTRACT 

HH 262 is a group of emitting knots displaying an "hour-glass" morphology in the Ha 
and [S 11] lines, located 3.5 arcmin to the northeast of the young stellar object L1551-IRS5, 
in Taurus. We present new results of the kinematics and physical conditions of HH 262 based 
on Integral Field Spectroscopy covering a field of ~ 1.5 x 3 arcmin^, which includes all the 
bright knots in HH 262. These data show complex kinematics and significant variations in 
physical conditions over the mapped region of HH 262 on a spatial scale of ^ 3 arcsec. A 
new result derived from the IFS data is the weakness of the [N ll] emission (below detection 
Hmit in most of the mapped region of HH 262), including the brightest central knots. Our 
data reinforce the association of HH 262 with the redshifted lobe of the evolved molecular 
outflow L1551-IRS5. The interaction of this outflow with a younger one, powered by L 1551 
NE, around the position of HH 262 could give rise to the complex morphology and kinematics 
ofHH 262. 

Key words: ISM: jets and outflows - ISM; individual: HH 262, L1551-IRS5, L1551 NE 



1 INTRODUCTION 



HH 262 is the designation of a group of emitting knots, sur- 
rounded by more diffuse emission, located ~ 3.5 arcmin to the 
northeast of the source L1551-IRS5 in Taurus. The two brightest 
HH 262 c ondensations (nam ed GH9 and GHIO) we r e first reported 
by Rodriguez et al.l ( Il989h and iGraham & HeveJ jl990l) . In the 
Ha+[N II] and [S II] narrow-band images, HH 262 shows a rather 
chaotic, "hour-glass"-shaped morphology that is quite different 
from the coUimated structures found in most of the H erbig-Haro 
(HH) jets jGamavich. Noriega-Crespo & Greenlll992llL6pez et al.l 
1 19951) . Apparently, the morphology of HH 262 seems to trace a 
wind-blown cavity from an embedded source located close to the 
GHIO knot. However, no point source suitable for powering the 
HH 262 emission has been currently reported at this location. Since 
its discovery, the origin of HH 262, as well as its relationship with 
other outfl ows in the L1551 star-fo rming region, has been matter of 
discussion. [Rodriguez et"aiHl989h first suggested that the GH9 and 
GHIO knots were associated with the redshifted lobe of the molec- 
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ular outflow powered by the source L1551-IRS5. They proposed 
that GH9 and GHIO may be the redshifted counterpart of a chain 
of HH bright knots (which includes the well-known HH objects 28 
and 29) extending from L1551-IRS5 to the southwest up to ~ 5 ar- 
cmin, and suggested that all these optical knots traced the walls of a 
wind cavity evacuated by the CO outflow powered by L1551-IRS5. 

Later on, Grah am & RubinI ('l992) obtained a long-slit spec- 
trum across the GH9 and GHIO knots. From the Ha line, they 
derived redshifted radial velocities covering a range extending 
over ~ 100 km s^^and measure peak velocities of +40 km s 
and -1-90 km s~^ and average velocities of -1-20 km s^^and 
+A5 km s~^in the GH9 and GHIO knots, respectively. These results 
reinforced the suggestion that GH9 and GHIO were associated with 
the redshifted lobe of the CO outflo w powered by L155 1-IRS5. The 
kinematics of HH 262, obtained bv'Lopez et al.l ( 1 19981) from proper 
motion measurements of the brightest knots, and from radial veloc- 
ities derived from Fabry-Perot observations in the [S II] A 6716 A 
line, is also consistent with the interpretation of the HH 262 ori- 
gin mentioned before. However, iDevine. ReiDurth&B"aiivl ( ll999h 
were uncertain whether the exciting source of HH 262 is L1551- 
IRS5 or L1551 NE, another younger, embedde d proto-binary 
source located ^ 150 arcsec northeast of IRS5 jEmerson et al.l 
ll984lRodrfguez. Anglada & Rag j| 199^ . L1551 NE also drives a 
molecular CO outflow, so that it is likely that the outflows pow- 
ered by these two sources are interacting. Further observations from 
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Figure 1. Diagram of the HH 262 region siiowing the location of HHs 
and YSOs. The contours correspond to the blueshifted and redshifted CO 
emission from the three molecular outflows identified in the region (two 
outflows in the NE-SW direction, associated with t he YSOs L1551 IRS 5 
and L1551 NE, and the east-west outflow: IMoriartv-Schieven et al.ll20o3) . 
The straight lines are drawn along the outflow axes. 



IStoiimirovic et alj ( l2006h showed evidence of blueshifted molec- 
ular gas arising close to LI 551 NE that extends to the north- 
eas t of the source and encompass es H H 262. As pointed out 
by iDevine, Rei purth & Baifyl ( Il999h and iMoriartv-Schieven et al.l 
( l2006h . HH 262 could be related to either the L1551-IRS5 or L1551 
NE outflows, its rather chaotic morphology originating through the 
interaction between the molecular outflows powered by these two 
sources (see diagram of Fig[T). Mo re recent sub-mm conti nuum ob- 
servations at 850 fim by Moriartv-Schieven et al.l | |2006|) detected 
emission around HH 262 that were interpreted as produced by 
warm dust entrained by the molecular outflows. HH 262 is thus 
a suitable place to search for optical signatures of the interactions 
among supersonic outflows. In addition, because of its morphology, 
HH 262 is better suited for performing Integral Field Spectroscopy 
(IPS) than long-slit observations when the kinematics and physical 
conditions of the emission need to be analysed. 

With the aim of advancing in the study of HH 262, we in- 
cluded this target within a program of Integral Field Spectroscopy 
of Herbig-Haro objects using the Potsdam Multi- Aperture Spec- 
trophotometer (PMAS) in the wide-field IPU mode PPAK. We 
present in this paper new results on the kinematics and physical 
conditions of HH 262 obtained from this IPS observing program. 
The paper is organized as follows. The observations and data reduc- 
tion are described in § 2. Results are given in § 3: The integrated 
maps in § 3.1, the channel maps in §3.2, and the spectral atlas of 
HH 262, including discussions of the characteristic of the more rel- 
evant knots in § 3.3. A global discussion and the main conclusions 
are given § 4. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations of HH 262 were made on 22 November 2004 with 
the 3.5-m telescope of the Calar Alto Observatory (CAHA). Data 
were acquired with the Integral Field Instrument Potsdam Multi- 
Aperture Spectrophotometer PMAS (Rothetal. 2005) using the 
PPAK configuration that has 33 1 science fibres, covering an hexag- 
onal FOV of 74x65 arcsec^ with a spatial sampling of 2.7 arcsec 
pe r fibre, and 36 a dditional fibres to sample the sky (see Fig. 5 
in lKelz et al.ll200d) . The 11200 grating was used, giving an effec- 



tive sampling of 0.3 A pix~^ (~ 15 km s^^ for Hq) and cover- 
ing the wavelength range ~ 6500-7000 A, thus including char- 
acteristic HH emission lines in this wavelength range (Hq, [N II] 
AA6548, 6584 A and [Sll] AA6716, 6731 A). The spectral reso- 
lution (i.e. instrumental profile) is ~ 2 A FWHM (~ 90 km s~^) 
and the accuracy in the determination of the position of the line cen- 
troid is ~ 0.2 A (~ 10 km s^^ for the strong observed emission 
lines). Eight overlapped pointings of 1800-s exposure time each 
were observed to obtain a mosaic of ~ 1.5 x 3 arcmin^ to cover 
almost the entire emission from HH 262. Note, however, that this 
mosaic does not cover the extended, arc-shaped emission known 
as HH 262E, which is most probably related to HH 262. In Ta- 
ble 1 we list the centre positions of each pointing and the cor- 
respon ding HH 262 knot identification following the nomencla- 
ture of [ Gamavich. Noriega-Crespo & GreenI ( 1 19921) . extended by 
iLopez etal.l ( ll995h . 

Data reducti on was pe rformed using a preliminary version of 
the R3D software (Sanchez 2006), in combination with IRAlQ and 
the Euro3D packages t Sanchez. 2004) . The reduction consists of the 
standard steps for fibre-based integral field spectroscopy. A mas- 
ter bias frame was created by averaging all the bias frames ob- 
served during the night and subtracted from the science frames. 
The location of the spectra in the CCD was determined using a 
continuum-illuminated exposure taken before the science expo- 
sures. Each spectrum was extracted from the science frames by 
co-adding the flux within an aperture of 5 pixels along the cross- 
dispersion axis for each pixel in th e dispersion ax is, and stored in 
a row-stacked-spectrum (RSS) file dSanchej^oioi) . At the date of 
the observations we lacked of a proper calibration unit for PPAK. 
Therefore the wavelength calibra tion was performe d iteratively, 
following all the steps described in iLopez et al.l ([2008') for HHL 73, 
another target of our HH IPS survey. The final wavelength solu- 
tion was set by using the sky emission lines found in the observed 
wavelength range. The accuracy achieved for the wavelength cali- 
bration was better than ~ 0.1 A (~ 5 km s^^). Observations of 
a standard star were used to perform a relative flux calibration. 
A final datacube containing the 2D spatial plus the spectral in- 
formation of HH 262 was then created from the 3D data by us- 
ing Euro3D tasks to interpolate the data spatially until reaching a 
final grid of 2.16 arcsec of spatial sampling and with a spectral 
sampling of 0.3 A . Further manipulation of this datacube, de- 
voted to obtain integrated emission line maps, channel maps and 
position-velocity maps, were made using several common-users 
tasks of STARLI NK, IRAF and GI LDAS astronomical packages and 
IDA ( Garcia-Lor enzo. Acosta-Pu lido & Megias- Fernandez 
an specific IDL software to analyse 3D data. 



3 RESULTS 

3.1 Integrated line-emission maps 

3.1.1 Emission line images 

Emission from Ha and [Sll] AA6716, 6731 A was detected in all 
the mapped region. In contrast, emission from the [N II] AA6548, 
6584 A lines was not detected, except in the knot GH9, towards 
the northwestern side of HH 262. Narrow-band images for each 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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Table 1. Positions of the single pointings of the HH 262 mosaic. 
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Figure 2. Integrated maps of the HH 262 emission obtained from the dat- 
acube by integrating the signal within the wavelength range including the 
hne labelled in each panel: (Ho : A6558-6569 A, [Nil] : A6580-6586 A, 
[S II] : A6716 6712-6719 A and [S II] : A6731 6728-6734 A). Fluxes have 
been displayed in a logarithmic scale, in units of 10~^^ erg s~^ cm~^. In 
all the maps, north is up and East is to the left. The HH 262 knots have been 
labelled in the Ha panel 



of the lines were generated from the datacube. For each position, 
the flux of the line was obtained by integrating the signal over the 
whole wavelength range covered by the line and subtracting the 
continuum emission obtained from the wavelengths adjacent to the 
line. The narrow-band maps of HH 262 obtained in this way are 
shown in Fig|2] 

As can be seen from the panels of the figure, the Hq and [S II] 
emissions reveal a similar morphology, although the emission of 
some the knots is stronger in [S II] than in Ha, while the brightness 
of the [S II] and Ha emission is similar for the other knots, in good 
agree ment with the results found in previous CCD narrow-band im- 
ages jGamavich. Noriega-Crespo & Greenlll992h . 

Our IPS data allow us to obtain the first [N II] narrow-band im- 





Figure 3. Velocity maps of the line centroid. Velocities (km are re- 
fen'ed to the Vlsr of the parent cloud. Contours of the integrated flux of 
the Ha line have been overlaid to help with the knots identification. 



age of HH 262. Because there are no CCD images acquired through 
a line filter with a passband narrow enough to isolate the [N II] lines 
from Hq, the morphology of HH 262 in the [N II] lines has re- 
mained unknown up to now. Note that the morphology of HH 262 
in [N II] is different than in Ha and [S II] : i. e. HH 262 does not 
show its characteristic "hour-glass" structure because the emission 
from [N II] was detected only in knot GH9. In the other knots, [N II] 
emission was only marginally detected or not detected at all. 



3.1.2 Kinematics 

Radial velocitjlf] maps for the Hq and [S II] emissions were ob- 
tained from the line centroids of a Gaussian fit to the line emis- 
sion. All the HH 262 knots show redshifted velocities in both lines 
(Hq and [S II]), the velocity field being rather complex (Figure|3]l. 
The velocities derived for the Ha emission range from ~ +25 to 
~ +85 km s~^ and, in general, are ^ 20 km lower than those 
derived from [S II] emission, which range from ~ +40 to ~ +110 
km s~^, the difference found between the Ha and [S II] velocities 
being reliable. The highest velocities are found towards the central 
knots of HH 262 (a velocity of ~ +70 in Hq and and ~ +100 
km s~^ in the [Sll] A 6716 A line, respectively, is found around 
GHIOW). The velocity decreases as we move north and south of 
the centre, with a steeper decrease towards the north. The lowest 
velocities appear towards the northwest, in the GH9 knots, where a 
value of ~ +30 km s^ ^ for both lines is derived around its northern 
edge. It is interesting to note here that the velocities derived in this 
work are in good agreement wit h those found in previo us studies 
(from long-slit spectroscopy of Graham & Rubinlll992l . and from 
Fabry-Perot observations of Lopez et al. 19981) . 

The FWHMs of the Ha and [S II] lines, obtained from the 
Gaussian fit and corrected from the instrumental profile, show the 
same behaviour: for both emissions, the widest FWHM values are 
found towards the centre of HH 262 (~ 80 km s~^ around the 
knot GHIOE and ~ 100 km s~^ around GHIOW), the line widths 
diminishing when we move towards the north, reaching values of 
~ 30 km s"^ around the knot 18. The FWHlVIs vary along GH9, 
ranging from ~ 40 km s^^ (similar to those found around knots 
19 and 20) in the southern region of GH9, and increase towards 



^ All the velocities in the paper are referred to the local standard of rest 
(LS R) frame. A Vlsr = +7 .0 km s~^ for the parent cloud has been taken 



from lOraham & RubinI h992h . 
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the north (a FWHM of ' 
edge). 



65 km s is found around the GH9 NE 



Table 2. Derived values for the [Nil] knots emission 



3.1.3 Physical conditions 

Line-ratio maps of the integrated line emission were derived to 
explore the spatial behaviour of the electron density (from [S II] 
6716/6731) and the excitation conditions (from [S Il]/Ha and from 
[N Il]/Ha at the loci where [N II] emission was detected). 

As a general trend, the electron density (ric) derived in HH 262 
is low when compared with densities found in most of the HH 
jets. The highest ric values are found in GH9, coinciding with the 
HH 262 loci where emission from [N II] was detected (ric = 490- 
1000 cm~^, increasing from the southeast to the northwest along 
GH9). The electron density decreases as we move towards the cen- 
tre of HH 262, from north (Wc = 270 cm"^ around K18) to south 
(ne = 100 cm"^ around K20). For the central knots (GHIOE, 
GHIOW and K25), the electron density lies close to the low-density 
limit of the [S II] 6716/6731 line ratio sensitivity and thus only an 
upper limit for ric of 30-50 cm~"^ was estimated. The electron den- 
sity increases slightly as we move from the centre towards the south 
(rie = 120cm^'^ around K22) and southeast (ric = 230 cm^'^ 
around K21). These values were derived from the [S II] 6716/6731 
ratio, using the TEMDEN task of the IRAF/STSDAS package, and 
assuming Tc — 10* K. 

The derived spatial distribution of the [S Il]/Ha line ratio indi- 
cates that the excitation increases from east to west across HH 262. 
The lower line ratio values ([S II]/Hq ~0.6-0.9) are found towards 
the western region of knots GHIOW and GH9, reaching the low- 
est value (highest excitation) around GH9, coinciding with the loci 
where [N II] emission is detected. Furthermore, the trend found is 
the [S II]/Hq line ratio decreases from the northern K 18-20 knots 
(with [Sll]/Ha ~ 1.8-2.2) to the centre of HH 262 (GHIOW-E 
knots, where [S Il]/Ha ~ 0.7-1.0). Beyond the centre of HH 262, 
the [S Il]/Ha line ratio increases slightly (the gas excitation de- 
creases) as we move to the southern knots (K22, where [S Il]/Ha ~ 
1.2). 

As already mentioned, the [N ll] emission appears very weak 
and it was detected only marginally in the brightest HH 262 cen- 
tral knots (see Fig.|2}. To provide better confirmation of this result, 
we have analysed an integrated spectrum of each knot, obtained by 
co-adding the signal of the spaxels inside the aperture encircling 
the corresponding knot emission (see Section 3.3), thus improving 
the SNR and the depth of the spectra. From these spectra, we have 
derived values for the [N II] flux emission (or an upper limit). Table 
2 lists the flux of the [N II] emission relative to the rms of the spec- 
trum, the [Nll]/Ha line ratio and its estimated error percentage. 
From the table, we can see that [N Il]/Ha line ratios of ~0.2 were 
obtained along GH9 knot. For the northern knots (K18, K19 and 
K20), we estimated an upper limit of ~0. 1 for this ratio (by consid- 
ering the [N II] flux emission enhancement at a 3a rms level). For 
the central (GHIO) and southern (K21, K26 and K22) knots, [Nil] 
flux emission over the 3a rms level was not detected at all. 

We have searched the literature looking for HH spectra that 
would be similar to those extracted for the HH 262 knots. Low 
[Nll]/Ha line ratios (~ 0.1-0.2), similar to those derived from 
the spectra of the HH 262 knots K19, K20 and GH9, have 
been measured in the spectra of kno wn HH jets (e. g. H H 49,50: 
Schwartz & Dopita 1980; HH 1 1 IL: iMorse et al.ll 19931: HH 123 : 
Reipurth & Heathcote, il99Q: H H 126B: lOgura & WalshI Il99ll : 
HH 34 (apex): Morse et aljl993h . However, the spectra of these jet 
knots, covering a wider wavelength range to include other emission 
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lines, seem to be quite different when compared with the HH 262 
extracted spectra, and would thus trace gas at different conditions. 
For example, the [S Il]/Ha ratios in these jet knots are in general 
lower than the [S II]/Hq ratio derived in the HH 262 knots having 
similar [N Il]/Ha values, and the derived of these jet knots are 
also higher than the densities derived for HH 262. Moreover, most 
of the above-mentioned knot spectra show emission from [O III] 
lines. Although the IFS data of HH 262 analysed here do not cover 
a wavelength range that includes the [O III] lines, it should be men- 
tioned that we failed to detect [Olll] emission at the location of 
HH 262 in a CCD image of 3000-s exposure time, acquired through 
an appropriate narrow-band filter in a previous observing run. Thus, 
the HH 262 spectra appear to have lower excitation than the knots 
mentioned before having similar [N II]/Hq: line ratios. 

In order to compare the line ratios obtained in HH 262 with 
current models of shock-ionization, it has to be noted that the 
HH 262 fluxes have not been corrected for reddening. Thus the 
derived line ratios are somewhat affected by differential extinc- 
tion. However, since the wavelengths of the lines considered are 
close, we expect that the differences between the observed and the 
corrected line ratios are not critical for model comparison. For in- 
stance, for other HH objects where corrected fluxes are available, 
the differences found are lower than 1 % and 15 % for the [N II]/Hq 
and [S Il]/Ha ratios respectively (see e. g. iRaga. Bohm & Cantdl 
1 19961 and references therein). 

Line ratios such as those derived from the HH 262 spectra can 
be reproduced in a framework of shocks with slow shock veloci- 
ties. To pr edict the expected line ratios of character istic HH emis- 
sion lines, iHartigan. Raymond & HartmanrJ l ll987l) modelled HH 
bow shocks for a wide range of shock velocities (14), pre-shock 
densities and pre-shock ionization conditions, assuming radiative 
planar shocks. Equilibrium preionization models in the low density 
regime (ric — 100-1000 cm^'') and with low shock velocities (14 
40 km s"^ ) give low [Nll]/Ha line ratios (~ 0.01-0.22, in- 
creasing with Vs, with a weak dependence on the electron density) 
and high [S Il]/Ha line ratios (~ 2.5-0.8, decreasing with Vs). Thus 
[N Il]/Ha (~ 0.2) and [S II]/Hq (~ 0.6) line ratios, such as those 
found in the GH9 knot, are compatible with those obtained in shock 
models with Vs — 30-40 km s^^ occurring in a low density and par- 
tially ionized medium. Shock models with higher shock velocities 
(Vs ^ 60 km s~^) predict [S Il]/Ha line ratios much lower (by a 
factor of 4), not compatible with the values derived in HH 262, even 
considering that the fluxes were uncorrected for reddening. Lower 
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[N Il]/Ha line ratios (~ 0.01) are found in the planar shock models 
with lower velocities (K ~ 20 km s~^). The line ratio values de- 
rived in the central knots GHIOE-W show a trend that seems to be 
compatible, within the errors, with that found in shock models of 
very low K- The [S II]/Hq line ratios derived in HH 262 are also 
well reproduced by the low- velocity (K ^ 40 km s~^) models of 
iHartigan. Morse & Ravmona ( 1 1994) that include the effect of the 
magnetic field, for low B values (~ 30 /iG, in accordance with the 
current estimates for the Taurus Molecular Cloud. ,Nakamura & Lil 
I2OO8I) . and low ionization fractions (a;e ^ 0.1). 



3.2 Channel maps 

Since the HH 262 emission spreads over a wide range of velocities, 
we will explore whether the complex, knotty HH 262 morphology 
shown in the narrow-band images is also found in all the channel 
maps or, on the contrary, there are changes in morphology depend- 
ing on the velocity. 



3.2.1 Spatial distribution of the emission 

A more detailed sampling of the gas kinematics was obtained by 
slicing the Ha and [S II] emission of the datacube into a set of ve- 
locity channels. Each slice was made for a constant wavelength 
bin of 0.3 A, corresponding to a velocity bin of ~ 14 km s^^ 
(13.7 km s^^ and 13.4 km s^^ for Ha and [S II] respectively). The 
channel maps obtained for Ha and [S II] A6716 A are shown in 
Figs . |4] and|5]respectively. 

From these maps, it can be seen that the emission spreads over 
velocity ranges of ^ 290 km s^^ and ~ 240 km s~^ wide for 
the Ha and [Sll] A6716 A emissions, respectively. The HH 262 
"hour-glass" structure is found for redshifted velocities ranging 
from ~ +5 to ~ +100 km s~^ in Ha, and from ~ +30 to 
~ +125 km s~^ in [S II] A6716 A . Furthermore, the channel maps 
also show that the northern (K18, K19 and K20) and northwestern 
(GH9) knots have emission at lower velocities as compared with 
the other HH 262 knots. These knots show emission at blueshifted 
velocities (up to ~ —50 km s~^) in contrast with the southern 
HH 262 knots (K22), whose emission just begins to appear at ve- 
locities close to the rest velocity of the ambient cloud. Note in par- 
ticular that the bright knot GH9 does not show significant emission 
for velocities higher than ~ +100 km s^^, while the emission from 
the central (GHIOE-W) and southern (K22) knots reach velocities 
up to ~-l-180 km s^^ and ~-l-150 km s~^ respectively. For GH9 
and the northern knots, the strongest emission corresponds to the 
velocity channels ranging from ~ +19 to ~ +45 km s~^ in Ha, 
and from ~ +30 to ~ +55 km s^^ in [S II] . In contrast, for the 
central GHIOE-W knots, the strongest emission appears at higher 
velocities (in channels ranging from ~ +60 to ^ +90 km s^^ in 
Ha, and from ~ +70 to ~ +110 km s~^ in [S II] ). Thus, we found 
signatures of deceleration of the optical emission as we move to the 
north, away from the exciting source of the molecular outflow. Fur- 
thermore, an increase of the density in the northwestern GH9 knot 
relative to the rest of the HH 262 could also be contributing to the 
deceleration of the outflow and thus derive lower velocities for the 
GH9 region. 
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Figure 4. Ha channel maps centred on the Vlsr radial velocities (km s~^) 
labelled in each panel. The spatial scale is given in one of the panels. 
Fluxes have been displayed in a logaiithmic scale, in units of 10^^® 



erg s ^ cm ^ . 



3.2.2 Density and excitation conditions 

The behaviour of the density and excitation as a function of veloc- 
ity was explored by generating the appropriate line-ratio channel 
maps. 

We did not find evidence of density variations as a function 
of velocity in HH 262, as the [S II] 6716/6731 line-ratio channel 
maps show a similar pattern for all the velocity bins. The [S Il]/Ha 
line-ratio channel maps are displayed in Fig. [6] The maps in Fig. |6] 
show a trend of an increasing [S Il]/Ha line ratio (i. e. a decreasing 
excitation) with increasing velocity. Regarding the spatial distribu- 
tion, the general trend found for each velocity channel is a decrease 
in the line ratios (i. e. an increase in the excitation) from east to 
west across HH 262. In fact, the behaviour found in the integrated 
line-ratio map, where the highest excitation (lowest line-ratio val- 
ues) is found along the western edge of knots GH9 and GHIOW, is 
also found for each velocity channel. 



3.3 The spectral atlas of HH 262 

A more careful inspection of the HH 262 IFS narrow-band im- 
ages (Fig. [2} shows the very complex morpholohy already found in 
the CCD images with higher spatial resolution: most of the named 
HH 262 knots appear to be quite extended and can be resolved into 
substructures (e. g. GH9). As already mentioned, IFS data are well 
suited to characterizing the spectral properties of such as complex 
morphology in a more efficient way that by using long-slit obser- 
vations. To take advantadge of this, we performed a more detailed 
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Figure 5. Same as Fig.O but for the [S II] 6716 A emission line. 



Figure 6. The [S II]/Hq line-ratio maps derived from the channel maps at 
the velocities labelled in the con'esponding panel. 



Table 3. Apertures defined for the extraction of the knot spectra. 



knot 



Aperture Spectra 
(arcsec) X (arcsec) 



K18 




15 Xl4 


46 


K19 




8.5x15 


27 


K20 




11 xl3 


30 




K20N 


4.5x6.5 


5 




K20S 


8.5x6.5 


12 


GH9 




13 X34.5 


96 




GH9N-shell 


2.5x24 


11 




GH9W 


5.5x7.5 


9 




GH9E 


8.5x15 


28 


GHIOW 




7.5x13 


20 


GHIOE 




17.5x14 


51 


K21 




7.5x4.5 


8 


K26 




9.5x5.5 


10 


K22 




11 xl3 


31 




K22N 


6.5x7.5 


11 




K22C 


6.5x8.5 


12 




K22S 


6.5x4.5 


5 



spectral atlas of HH 262. First of all, we obtained a representative 
spectrum for each knot by co-adding the IFS spectra extracted from 
the spaxels inside an appropriate aperture encircling the emission 
of the HH 262 knot. The size of the apertures and the number of 
spectra that were co-added to obtain the representative spectrum of 
the knot, are listed in Table 3. 



Fig. |7] displays the spectra of the HH 262 knots, obtained in 
such a way, in the spectral windows covering the Ha + [N II] (left 
panel) and [S II] (right panel) wavelength ranges. In order to make 
easier comparisons between the lines in the two spectral windows, 
the intensities in each of the spectra were normalized to the corre- 
sponding Ha peak intensity. A quick inspection of the spectra in the 
two spectral windows suggests that the excitation conditions signif- 
icantly vary from knot to knot. Note in particular that the spectrum 
extracted at the northern knot (K18) presents the highest [S Il]/Ha 
line ratio (this indicating the lowest gas excitation). Concerning the 
line profiles, a more careful inspection of the spectra suggests that 
the line profiles also vary from knot to knot. In particular, it can 
be appreciated that the southern (K21, K22 and K26) and central 
(GHIOW-E) knots show more complex, asymmetric line profiles 
than the northern (K18, K19) knots for both the Ha and [S II] lines. 

Furthermore, in order to better visualize the spectral variations 
at a smaller spatial scale, both along and across HH 262, we ex- 
tracted spectra of individual spaxels for two cuts made along the 
two main HH 262 axes. One cut was performed from north to 
south, intersecting emission from knots K18, KI9, K20, GHIOE 
and K22. The other cut, from northwest to southeast, intersects 
emission from knots GH9, GHIOW, GHIOE and K21. The position 
where the Ha emission has the peak, at the central knot GHIOE, 
was set as the origin of positions for referring the location of a given 
spectrum. The resulting spectra from these two cuts, in the [S II] 
window and normalized to their Ha peak intensities, are displayed 
in Figs.[8]and|9] The spectra displayed in these figures suggest there 
are variations in the kinematics and physical conditions (excitation 
and density) along and across HH 262 even on a short spatial scale 



Integral Field Spectroscopy ofHH 262 7 




Figure 7. Spectra of the HH 262 knots in the spectral windows including Ha+[NlI] (a) and [SlI] 6716,6731 emission lines (c). Each spectrum has been 
obtained by co-adding the spectra of the spaxels inside an aperture (continuous box) including the knot emission, and its intensity has been normalized to 
the Ho peak intensity. The apertures have been marked and labelled over the Ho map of HH 262 (b) and listed in Table 3. Dashed boxes correspond to the 
substructures defined inside a given knot. 



3 arcsec), because we found appreciable differences, both in the 
hne intensity and in the line profiles between the spectra extracted 
at adjacent spaxels. In particular, the longitudinal cut (Fig.[8]l shows 
the already mentioned increasing in the gas excitation as we move 
from the north to the centre of HH 262, suggested by the decrease 
along this direction in the [S II] line intensities relative to Ha. This 
cut also shows that, as a general trend, the line profiles for positions 
north of the central HH 262 knots are more symmetric and single- 
peaked as compared with those of the centre and south, where the 
profiles appear wider, suggesting in some cases a contribution from 
two velocity components, partially resolved at our spectral resolu- 
tion. Such a regular trend is not observed in the transversal cut (Fig. 
|9]l, although the extracted spectra reinforced the finding of the short 
spatial-scale variation conditions across HH 262. 

3.3.1 Remarks on selected knots 

As has already been mentioned, the Ha and [S II] narrow-band im- 
ages of HH 262 show that the emission from some knots appears 
split into subcondensations. Furthermore, some knots (e. g. K22 
and the central GHIOE-W knots) show more complex, asymmetric 
line profiles, suggesting a contribution from two velocity compo- 
nents at several positions, while the line profiles are found to be 
more symmetric and single-peaked in the other knots. Thus, in or- 
der to characterize better the kinematics and physical conditions 
in selected HH 262 knots, we also extracted spectra for each knot 
subcondensation (see Table 3 and Fig. [T] central panel for identi- 
fication) in addition to the spectrum extracted for the whole knot 
emission. 

• K20: 

Two subcondensations were morphologically defined in knot 20 
(K20-N and K20-S), its peaks being found to be separated by ~ 
10 arcsec. We have not found appreciable kinematic differences 
between these two subcondensations. As can be seen in Fig. [TO] 
the Ha and [S II] line profiles are both similar for the two subcon- 
densations as well as the velocities derived from the centroid of 
these lines (^^-1-60 km s~^ and -1-85 km s~^ respectively). In con- 



trast, a significant change in the [S II]/Hq line ratios can be appre- 
ciated when comparing the spectra of the two K20 subcondensa- 
tions, which traces an appreciable increase in the gas excitation as 
we move from north to south by only a few arcsec, suggesting a 
steep gradient of the gas excitation along K20. 

• GH9: 

The IFS data indicate that the northwestern arm of HH 262 (i. e. 
GH9) seems to have different physical conditions compared to the 
other HH 262 knots. We extracted the spectra for two subcon- 
densations (GH9W and GH9E) and also for the lower-brightness 
shell at the northern GH9 edge (N-Shell) (see Fig. |7] and Table 3 
for identification). Their spectra are shown in Fig. [12] Concern- 
ing the kinematics, we did not find appreciable differences among 
the GH9 subcondensations. The line profiles appear to be single- 
peaked along the whole of GH9. The velocities derived from the 
line centroids are similar in the N-Shell and GH9W, being slightly 
higher (by 15 km s^^) in GH9E. In contrast, we can appre- 
ciate changes in the physical conditions within GH9 as we move 
along the knot. Differences in the excitation conditions (increas- 
ing from N-Shell, GH9E to GH9W) are present as the [S II]/Hq 
line ratio variations indicate. Furthermore, the electron density also 
varies from one subcondensation to another. The highest rio (~ 
1000 cm~^) is found in the N-Shell, decreasing to ~ 850 cm"^) 
in GH9W and to ric ~ 600 cm~^ in GH9E. Taking into account that 
71c drop at the centre of HH 262 (i. e. around knots GHIOE-W) until 
reaching values for which the [S II] line ratio is not density sensi- 
tive in the central region, we can infer that there is a steep density 
gradient across HH 262 in the northwest direction. 

• K22: 

K22 is one of the HH 262 knots whose extracted spectrum sug- 
gests that the line profiles are double-peaked. Morphologically, the 
knot appears curved, elongated and quite extended, but split into 
subcondensations. We extracted spectra for each of the three bright 
K22 subcondensations (K22-N, K22-C and K22-S, the separation 
between the K22-N and K22-S peaks being ~ 10 arcsec). The spec- 
tra are shown in Fig. [TT] There are appreciable differences in the 
[S Il]/Ha line ratios among the K22 subcondensations, this being 
indicative of an increase in the gas excitation as we move from 
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Figure 8. HH 262 spectra in the [S II] spectral window extracted from in- 
dividual spaxels along a longitudinal axis. Spectra are displayed from north 
(top left) to south (bottom right). Each spectrum corresponds to an aperture 
of 2.16 X 2.16 arcsec^ and has been normalized to the Ha peak intensity. 
The knot where the spectrum was extracted is indicated in the correspond- 
ing panel. For each spectrum, the position within HH 262 is indicated by 
the (a, S) offset in pixel units (2.16 arcsec pixel^^) from the reference po- 
sition. The reference (0,0) position has been set in the Ho peak position of 
GHIOE. 



south to north along K22. Thus, a similar gradient in the gas exci- 
tation, but in the opposite direction to that found in K20, is detected 
along K22. The rtc derived from the [S II] line ratio increases from 
south to north along K22, i. e. following the same trend that the gas 
excitation. 

Concerning the kinematics, both the Ha and [S II] line profiles 
of the three subcondensations appear asymmetric and would appear 
to be double-peaked. We used the [S II] lines to perform a decom- 
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. Same as Fig.[8] but for the spaxels of a cut along the transversal, 
southeast HH 262 axis (from top left to bottom right). 
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Figure 10. Spectra in the Ha+[N 11] (left) and [S II] (right) windows of 
the two K20 knot subcondensations obtained by co-adding the IPS spectra 
extracted from the spaxels included in the apertures encircling the subcon- 
densations. In the two windows, the intensity has been normalized to the 
con'esponding Ha peak intensity. 



position of the line profiles by applying a model with two Gaus- 
sian components. Gaussian fitting was perfomed using the DIPSO 
package of the STARLINK astronomical software. The wavelength 
difference between the two [S II] lines was fixed according to their 
corresponding atomic parameters. The same FWHM of the Gaus- 
sian in each of the two [S II] lines was assumed for each component 
of the profile. According to the fit parameters, we found two veloc- 
ity components, both redshifted and with a peak separation of ~ 
50 km s~^ between them. For both components, the trend found is 
the peak velocity increases from the southern to the northern sub- 
condensation. The strongest component has lower redshifted peak 
velocities (Vlsr from -1-55 to -1-70 km s~^ from south to north), 
while the faintest component has higher redshifted peak velocities 



Integral Field Spectroscopy ofHH 262 9 



CO 

CL, 

Q 



0) 




- 1 



6560 



6580 



X (A) 



6720 

A (A) 



6740 



Figure 11. Same as FigfTolbut for K22. The fits to the [S II] line profiles, ob- 
tained from the model, have been drawn in the right panel. The two velocity 
components are shown at the bottom. The residual, obtained by substracting 
the fit from the observed profile, is plotted in green. 





' 1 


, , 1 

GH9-shell -n 




1 1 1 1 1 

GH9-sheIl - 


( Peak 










o 




H 1 1 1 \ 


1 1 1 1 1 1 1 







GH9-W - 




GH9-W - 


>> 










itensi 








\ A 


)-< 






.ve 


1 


1 GH9-E - 


1 1 


1 1 1 1 1 

GH9-E - 


Relati 



















6560 6580 

X (A) 



6720 6740 

X (A) 



Figure 12. Same as Fig.[To]but for GH9 knot. 

(Vlsr from +105 to +120 km from south to north). In Fig. 
[TT] (right panel), we have also drawn the fits to the [S II] line pro- 
files obtained from the model. The two individual components are 
shown at the bottom. Residuals, obtained by subtracting the fit from 
the observed line, are also plotted. 
• GHIO: 

Fig. [Ts] displays the spectra of the central GHIOE (top) and 
GHIOW (bottom) knots in the Ha+[N II] (left) and [S II] windows. 
In the entire area analysed, the spectra extracted at each spaxel of 
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Figure 13. Same as Fig.[TT]but for the GHIOE-W knots. 



the region covering these two central knots show broad and asym- 
metric line profiles, suggesting in most of the positions a double- 
peaked profile. We performed the line profile decomposition analy- 
sis described above (K22 subsection) in all the spaxels of this cen- 
tral region using the [S II] doublet. From the fit, we found two red- 
shifted velocity components spreading over the whole region, the 
component with the lower velocity being the strongest. We were 
unable to find a clear spatial separation between these two com- 
ponents that would allow the association of each of them with a 
given region of the central knots. However, a certain trend appears 
in which the higher-redshifted component is brighter towards the 
east of the GHIOE knot, having its peak intensity somewhat dis- 
placed from the knot centre and showing in addition a secondary 
peak intensity in a region around the centre of GHIOW. The lower- 
redshifted velocity component has the peak intensity in the regions 
located around the centre and southeast of GHIOE. 

We show in Fig. [T3] the fits to the [S II] line profiles obtained 
from the model applied to the spectra, obtained as described in §3.3, 
of knots GHIOE and GHIOW. From the fit parameters, we found 
Vlsr of +60 and +120 km s^^ for the two GHIOE components 
and of +85 and +145 km s~^ for these of GHIOW knot. 



4 DISCUSSION AND CONCLUSIONS 

We have carried out a study of the kinematics and physical condi- 
tions of HH 262 based on IFS observations. From the analysis of 
these data we find: 

• HH 262 shows a non-coUimated, "hour-glass" morphology in 
the Ha and [S II] lines. In contrast, the emission from the [N II] 
lines is below detection limit except in the northwestern knot GH9. 
This work therefore presents the first narrow-band [N II] image cur- 
rently obtained of HH 262 to show its surprisingly different mor- 
phology as compared with the morpholgy found in the other char- 
acteristic HH emission lines. 

• In all the HH 262 knots, radial velocities derived from the line 
centroids appear redshifted, relative to the cloud velocity. The high- 
est velocities and FWHM values are found in the central knots. 

• A more detailed picture of the kinematics was derived from 
the channel maps, which shows that the emission spreads over a 
wide range of velocities, the range slightly varying for the differ- 
ent knots. Furthermore, inspection of the HH 262 extracted spectra 
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shows that there are broad, asymmetric hne profiles at the south- 
ern (K22) and central knots (GHIOE-W) that are double-peaked at 
several positions, which would suggest contribution from two ve- 
locity components. The line profiles appear more symmetric and 
single-peaked away from the centre, in the northern (K18-K20) 
and northwestern (GH9) knots. 

• Concerning the physical conditions in HH 262, we found vari- 
ations on a spatial scale of ^3 arcsec (i. e. of the order of our spa- 
tial resolution). As a general trend, the electron densities, derived 
from the [S II] line ratios, are low. In particular, only an upper limit 
could be set for the brightest, central HH 262 knots. Because ric is 
referred to the density of the ionized gas, we propose that the only a 
small fraction of the atomic gas in HH 262 is ionized, the remainder 
being mostly neutral gas. It is also plausible that the total density 
of the gas is lower in this region because the evolved molecular 
outflow has had enough time to clear the medium. Over the entire 
velocity range of the emission, the gas excitation, derived from the 
[S Il]/Ha line ratios, increases as we move towards the western re- 
gion of knot GHIOW and to knot GH9 i.e. moving towards the 
edges of HH 262 which are closer to the higher extinction region 
of the cloud. The general trend found is thus that the densest knots 
have the highest gas excitation and the lower centroid velocities. 

All these characteristics of HH 262 derived from the IPS data 
reinforce the causal relationship of the optical line emission with 
the two molecular outflows (one powered by L1551-IRS5 and the 
other, by L1551 NE) found in the region were HH 262 is located. 
In fact, is widely accepted the association between HH 262 and 
the redshifted lobe of the L155I-I RS5 molecular outflow, first pro- 
posed by Rodrigue z et al I h989l) and confi rmed later on by the 
Fabry-Perot observations 
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Detailed P V maps of the L1551-IRS5 outflow performed 
by IStoiimirovic et al.l j2006l) from high-sensitivity CO observa- 
tions show a pattern that is characteristic of an evolved molecu- 
lar outflow, according to c urrent model simulations of the entrain- 
ment outflow mechanisms jVelusamv & Langedl 19981 .1 Richer et al.l 
l200d) . In the lStoiimirovic et alj j2006h PV CO maps, the HH ob- 
jects are found associated with the molecular emission at the high- 
est velocities, highlighting the interaction regions of the shocks at 
the head of the shock interface between the supersonic gas and 
the environment. In contrast, the CO emission at lower velocities 
appears as a lim b-brightened shell-like em ission. Thus, the result 
derived from the IStoiimirovic elTa work support the pro- 

posed scenario, in which the HH 262 emission traces the walls of a 
conical cavity evacuated by an evolved molecular outflow (L1551- 
IRS5). In such a scenario, one would expect that the HH 262 radial 
velocities to decrease (due to projection effects) as we move away 
from the axis outflow (i. e. from the central GHIO knots, where we 
found the highest integrated velocities, to the other knots). It is also 
expected to find double-peaked line profiles along the outflow axis 
(corresponding to material moving in the forward and rear faces of 
the cone) and a single-peaked profile away from the axis. Such a 
trend is observed in HH 262 from our IPS data, where the lines of 
the central knots appear with a complex, asymmetric line profiles, 
being double-peaked for most positions, while the line profiles of 
the knots farther away from the outflow axis (e.g. K18, K19 or 
GH9) appear more symmetric and single-peaked. 

In addition to the L1551-IRS5 outflow, the blueshifted lobe 
of another younger, dimme r CO outflow, powered by L15 51 NE, 
also encompases HH 262 iMoriartv-Schieven et al ] i2006h . Thus, 
HH 262 is located in a region where two molecular outflows are 
interacting (see Fig[T) and it is expected that this interaction gives 



signatures on the jet turbulence. This interaction could in part be 
responsible for the complex velocity pattern and the wide line pro- 
files found in the optical emission lines. In a turbulent jet scenario, 
the highest PWHM values are expected to be found close to the 
axis, where there are contributions from gas within a wider range 
of velocities because we intersect the entire jet beam, while PWHM 
values would be lower at the edge of the jet, away from the axis. 
This is the trend found in the IPS HH 262 data for the PWHMs , 
where the highest values correspond to the central GHIO knots and 
decrease moving towards the edges of the emission. 

The association of HH 262 with an evolved molecular outflow 
near the edge of a dark cloud would in part justify the low density 
values derived for most of the knots, as the outflow has cleared 
the medium. Note however that the density of the gas is probably 
higher than the density derived from the [S II] line ratio, which is 
an indicator of only the density of the ionized fraction of the gas 
rather than of the total density ( Sanchez et al. 2007). 

Pinally, let us discuss the new finding presented here con- 
cerning the low intensity of the [N II] emission in HH 262. Such 
a behaviour has been found in slow-velocity shoc ks as those 
modelled bv 'H artigan. Raymond & HartmanrJ ( Il987l) . In particu- 
lar, the [N Il]/Ha and [S II]/Hq line ratios derived for the knot 
GH9 are fully compatible with those derived in shock models 
with Vs ~ 30^0 km s~^ occurring in a low density, partially 
ionized medium. Lower [Nll]/Ha line ratios (~ O.OI) are found 
in models with lower Vb (^ 20 km s~^ ), thus being compati- 
ble, within the errors, with the marginal detection of [N II] in the 
northern knots and with the line ratios derived in those knots. 
The HH 262 line ratios are also compatible with thos e predi cted 
by the shock models of iHartigan. Morse & Ravmonj ( 1 19941) for 
low magnetic field strength and low ionization fraction values. 
Note however that the morphology of HH 262 is too complex 
and rather different from the morphologies found in most known 
HH jets. Thus, it should not be exp ected that the radiative pla- 
nar shock models of [H artigan, Ravm ond & Hartmamj ( Il987t) and 
IHartigan. Morse & Ravmond ( II994I) succeeded in reproducing all 
the HH 262 properties derived from the current IPS data. Purther 
IPS observations covering a wavelength range blueward of [N II] -l- 
Ha are needed to search for emission in the [Nl] A 5200 A line 
whose strength could justify the weakness of the [N II] emission 
found in HH 262 provided the nitrogen abundance were not far 
from the solar value. 
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